A NEW APPROACH TO THE 2- VARIABLE 
SUBNORMAL COMPLETION PROBLEM 
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Abstract. We study the Subnormal Completion Problem (SCP) for 2-variable weighted shifts. We use 
tools and techniques from the theory of truncated moment problems to give a general strategy to solve 
SCP. We then show that when all quadratic moments are known (equivalently, when the initial segment 
of weights consists of five independent data points), the natural necessary conditions for the existence of 
a subnormal completion are also sufficient. To calculate explicitly the associated Berger measure, we 
compute the algebraic variety of the associated truncated moment problem; it turns out that this algebraic 
variety is precisely the support of the Berger measure of the subnormal completion. 



1. Introduction 

We present a new approach to the Subnormal Completion Problem (SCP) for 2-variable weighted 
shifts. It employs the localizing matrices introduced and studied in [11] in the context of the truncated 
if-moment problem (if-TMP). This helps identify potential candidates for weights, and makes the 
problem more accessible. 

We first give a general strategy to solve SCP, and we later apply it to solve the SCP with quadratic 
data. That is, given an initial set of weights Q± consisting of five independent data points (aoo> Poo, &io, 
aoi and An); we prove that the natural necessary condition for the existence of a subnormal completion 
is also sufficient. Concretely, associated to the five given weights is a 3 x 3 moment matrix M(Q\), 
whose positive semi-definiteness is a necessary condition for the existence of a subnormal completion; 
in symbols, M(fti) := (7u+v) U) vez2 ,|u|,|v|<i3 wnere 7oo := 1, 7io := «oo> 7oi := Poo, 720 := «i «00' 
7li := Ooi/^OO' an d 702 := /^oi/^oo- We prove that the necessary condition M(Q\) > turns out to 
be sufficient for the existence of a representing measure \x supported in and satisfying the property 
supp [i n (0, +oo) 2 ^ 0; the measure fi then gives rise to a subnormal completion of fix. Once we know 
that a representing measure exists, we use techniques from the theory of truncated moment problems to 
find a concrete expression for it. 

As a first step, we build new weights 020, an, «02 and P02, an d we use them to construct the localizing 
matrices M X (Q^) and M y (^), where is a proposed extension of fix. The positive semi-definiteness of 
M(f2x) is then used to establish that the localizing matrices M x (Cl^) and M y {Vt^) can be made positive 
semi-definite for suitable choices of the new weights «20j a \i, a 02 and /3o2- That is, the condition 
M(fix) > triggers the two conditions M^f^) > and M y {Cl^) > for appropriate values of 0:20, a \i, 
CK02 and /3o2- Once that happens, we prove that a flat (i.e., rank-preserving) extension M(Q^) of M(fix) 
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exists, thereby giving rise to a unique representing measure [i for O3, which is the Berger measure of the 
subnormal completion. The explicit form of fi can be obtained by first determining the support of fj,, 
which agrees with the algebraic variety of ^3. 

In one variable, SCP was stated and solved in [7J: 

Problem 1.1. (One-Variable Subnormal Completion Problem) Given m > and a finite collection of 
positive numbers U m = {afc}£l ' fi n< ^ necessary and sufficient conditions on O m to guarantee the existence 
of a subnormal weighted shift whose initial weights are given by Q m . 

Since subnormality implies hyponormality, the condition ao < ol\ < a% < • • • < a m is obviously 
necessary; moreover, it is easy to dispose of the case when = a^+i for some < k < m — 1, so one 
can always assume that ao < o-i < ■ ■ ■ < a m . 

The cases m = and m = 1 are straightforward, with canonical completions given by ao, ao, ao, ■ ■ ■ 
and ao, a±, a%, ■ ■ ■ , respectively. The solution of the case m = 2 is based on the positivity of the moment 

and of the localizing matrix H r (2) := ; the explicit calculation 

7i 72 J 6 V ; V 72 73 /' y 

of the subnormal completion requires recursively generated weighted shifts Example 3.12]. 

In the general (1-variable) case, the Subnormal Completion Criterion (SCC) [7 S Theorem 3.5] states 
that a subnormal completion exists if and only if an £-hyponormal completion exists, where i := [^] + 1> 

Theorem 1.2. (One-Variable Subnormal Completion Criterion; cf. [3 Theorem 3.5] ) Let O m = {afclfcLo 
be a finite collection of positive numbers, let k := [^5—] and I := [y] + 1, and let H(k) = H(p, m ) := 
(li+j)o<i,j<k, H x (£-1) = H x {Vt m ) := (7j +i+ i)o<i i j<<?_i and := (7, 7 i+1 ••• H+j) 7 '■ The 

following statements are equivalent. 

(i) £l m admits a subnormal completion; 

(ii) £l m admits an t-hyponormal completion; 

(Hi) H(k) > 0, H x (£-1) > 0, and v(k + l,k) £ Ran H(k) if m is even (v(£ + l,£-l) £ Ran H x (£-1) 
if m is odd); 

(iv) H(Q m ) admits a positive flat (i.e., rank-preserving) extension H(Cl m+ i) such that H x (Cl m+ i) > 0. 

We now formulate the 2- variable SCP: 

Problem 1.3. (2-variable Subnormal Completion Problem) Given m > and a finite collection of pairs 
of positive numbers VL m = {(a^, /?k)}|k|<m satisfying \2. 1\) for all |k| < m (where |k| := k\ + k2), find 
necessary and sufficient conditions to guarantee the existence of a subnormal 2-variable weighted shift 
whose initial weights are given by £l m . 

While the research in [7] provided a complete solution to SCP in one variable, the 2-variable version 
requires new tools and techniques. At present, no general solution exists, and the problem appears to be 
quite difficult. When m = 0, in one variable the canonical subnormal completion of ao is the weighted 
shift ao,ao,ao, ■ ■ ., with Berger measure fi := 5 a 2; in two variables, the canonical subnormal completion 
of (ojoOj A)o) is the 2-variable weighted with weight sequences a^ := ctoo and := /3oo (all i, j > 0) and 

Berger measure a := 5 n 2 x 5 R 2 . 

r a oo Poo 

When m = 1, the 1-variable case is still straightforward; i.e., the canonical subnormal completion is 

2 2 

ao,a\,a\, . . ., with Berger measure (1 |)<5o H — ^S a 2. In two variables, however, the problem becomes 

highly nontrivial. For the singular case, and using the results in J9j Section 6], C. Li gave in [21] a solution, 
which seems a bit ad hoc and unmotivated, with extensive calculations using Mathematica. The proof in 
[9l pages 39 and 40] establishes the existence of a representing measure [i for SCP with quadratic moment 
data (this is the case m = 1 in two variables); however, the ensuing statement that supp /i C is made 
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without a proof, and it does not appear to follow easily from the comments preceding it. It is indeed 
true, as we show in the present paper using the tools and techniques from 

In Section [5] below, we shall apply our general strategy to solve SCP to the case m = 1 and prove that 
a representing measure always exist if the associated moment matrix M(l) is positive semi-definite. In 
Section [6] we shall calculate the Berger measure using canonical column relations in the flat extension 
M(2) of M(l). The reader will note how effective the theory of truncated moment problems can be in 
detecting the location of the atoms of the unique representing measure for M(2); this is in sharp contrast 
with the ad hoc techniques and extensive symbolic manipulation present in [21j. 

The case m = 2, in full generality, remains open; however, in Example 14 . 4 1 b elow we solve SCP whenever 
the associated moment matrix M(l) is singular. For m > 3, the results in [7] and [8] show that, in the 
1-variable case, it is not always possible to build a subnormal completion; of course the same is true in 
two variables: indeed, if ao,a±, 02,03 is a collection of weights admitting no subnormal completion, it 
suffices to consider the 2- variable collection given by := and /?k := 1 (|k| < 3) in order to produce 
such an example. 



Problem 11.31 is closely related to truncated moment problems. Given real numbers 7 



7 



(2n) 



700,701,710,702 , 7n,720, ■ ■ ■ ,702n, • • • ,72n0 with 700 > 0, the truncated real moment problem for 7 entails 
finding conditions for the existence of a positive Borel measure /i, supported in R 2 , such that 



y l x J dfi, < i + j < n. 



Given 7 = r y^ 2n \ we can build an associated moment matrix M(n) = M(n)(j) := {M[i,j]{ n f))2j—Q, where 



M[i,j]( 7 ) := 



/ 7o,i+j 7i,i+j-i 

7l,i+i-l 72,i+j-2 



V 7i,j 7i+lJ-l • 

We denote the successive rows and columns of M(n)( 7 ) by 



1,X,Y,X 2 ,YX, Y 2 , 



7j+l,i-l 

7i+j,o J 
,Y n . 



Observe that each block M[i,j](j) is of Hankel form, i.e., constant in cross-diagonals. (For basic results 
about truncated moment problems we refer to [9] and [11], ) 

We conclude this section by stating a result from |14j , which we will need in Section [3l Recall that a 
commuting pair {T\,T2) is 2-hyponormal if the 5-tuple (Tj, T2, T 2 , T\T%, T| ) is hyponormal (cf. Section 
[2] below). For 2- variable weighted shifts, this is equivalent to the condition 

l\) (cf. PH Theorem 2.4]), 



M u (2) := (j u+ (rn,n)+( P ,q))o<m+n<2 > (all U € 

0<p+q<2 



that is, 



/ 7u 7u+(0,l) 7u+(l,0) 7u+(0,2) 7u+(l,l) 7u+(2,0) \ 

7u+(0,l) 7u+(0,2) 7u+(l,l) 7u+(0,3) 7u+(l,2) 7u+(2,l) 

7u+(l,0) 7u+(l,l) 7u+(2,0) 7u+(l,2) 7u+(2,l) 7u+(3,0) 

7u+(0,2) 7u+(0,3) 7u+(l,2) 7u+(0,4) 7u+(l,3) 7u+(2,2) 

7u+(l,l) 7u+(l,2) 7u+(2,l) 7u+(l,3) 7u+(2,2) 7u+(3,l) 



> (all u € 



1.1) 



\ 7u+(2,0) 7u+(2,l) 7u+(3,0) 7u+(2,2) 7u+(3,l) 7u+(4,0) / 

An entirely similar formulation exists for ^-hyponormality (£ > 1), i.e., one requires 

M u {£) := ( 7 u+( m ,n)+( P , g ))o< m+ n<. > (all u G Z 2 ) (cf. [H Theorem 2.4]). 

0<p-\-q<£ 



1.2) 



2. Notation and Preliminaries 

Let Ti be a complex Hilbert space and let B(H) denote the algebra of bounded linear operators on 
rl. We say that T £ B{%) is normal if T*T = TT* , subnormal if T = N\u, where N is normal and 
N(H)Q H, and hyponormal if T*T > TT*. For S,T £ B(ft) let [S,T] := ST - TS. We say that an 
n-tuple T = (Ti, • • • ,T n ) of operators on T~L is (jointly) hyponormal if the operator matrix 



[T*,T1 



/ [T-f, Ti] [T|, Ti] ••• [T*,Ti] \ 

[2?,r 2 ] [t*,t 2 ] ••• [t*,t 2 ] 
V [rf,r n ] [t 2 *,t„] ••• [T*,T n ] J 



is positive semi-definite on the direct sum of n copies of Ti (cf. pQ, [E]). The n-tuple T is said to be 
normal if T is commuting and each Ti is normal, and T is subnormal if T is the restriction of a normal 
n-tuple to a common invariant subspace. Clearly, normal => subnormal => hyponormal. 

The Bram-Halmos criterion for sub normality states that an operator T £ B(rl) is subnormal if and 
only if 

^(T^-,T^)>0 

for all finite collections Xo,x±,--- ,Xk £ Ti ([3], [4]). Using Choleski's algorithm for operator matrices 
[22] . it is easy to see that this is equivalent to asserting that the fc-tuple (T, T 2 , • • • ,T fc ) is hyponormal 
for all k > 1. 

For k > 1, we say that a commuting pair T = (Ti,T2) is k-hyponormal if T(fc) := (Ti, T 2 , T 2 , T2T1, 
Tf, ••• ,Tf,T 2 T?-\--- ,T|) is hyponormal (cf. [13] )• Clearly, subnormal =^ (A; + l)-hyponormal 
/c-hyponormal for every k > 1, and of course 1-hyponormality agrees with the usual definition of joint 
hyponormality. The multivariable Bram-Halmos criterion was obtained in [T5j, and its formulation is 
essentially identical to the 1-variable one: T is subnormal if and only if T(/c) is hyponormal for all k > 1. 

For a = {a n }^L a bounded sequence of positive real numbers (called weights), let W a : £ 2 (Z + ) — > 
£ 2 (Z + ) be the associated unilateral weighted shift, defined by W a e n := a n e n+ \ (all n > 0), where {e n }^L 
is the canonical orthonormal basis in £ 2 (Z + ). The moments of a are given as 



7fc = lk{a) 



1 if fe = 

a?, ■ ■ ■ al , if A: > 



It is easy to see that VF a is never normal, and that it is hyponormal if and only if an < a\ < 
Similarly, consider double-indexed positive bounded sequences a = {ak},/3 = {/?k} £ ^^(Z 2 ^), k = 
(fci, fc 2 ) £ ^+ := x Z+ and let l 2 (Z\) be the Hilbert space of square-summable complex sequences 
indexed by (Recall that £ 2 (Z^_) is canonically isometrically isomorphic to £ 2 (Z + ) £ 2 (Z + ).) We 

define the 2- variable weighted shift T = (Ti,T2) by 



Tie k := a k e k+£l 
T 2 e k := /3 k e k+£2 , 

where e\ := (1,0) and e 2 := (0, 1). Clearly, 

T X T 2 = T 2 Ti ^ /3 k+£1 Q k = a k+£2 /3 k (all k £Z 2 ). (2.1) 



In an entirely similar way one can define multivariable weighted shifts. 
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Given k E Z+, the moment of (a, j3) of order k is 

1 if k = 



7k = 7k(a,/3) :-- 



a2 (o,o) ■ ■ ■ a (fc!-i,o) if fci > 1 and k 2 = 



£(0,0) ' ' ' /W 2 -i) if fci = and k 2 > 1 

a (o,o) • • • a (fex-i,o) • P(ki,°) ' ' ' ^(*i,*2-i) if kl - 1 and * 2 - 1 



We remark that, due to the commutativity condition (12, ip . 7^ can be computed using any nondecreasing 
path from (0,0) to (ki,k 2 )- 

We also recall a well known characterization of subnormality for multivariable weighted shifts [20J , due 
to C. Berger (and independently to R. Gellar and L.J. Wallen [19]) in the 1-variable case: T = (Ti, • • • , T n ) 
is subnormal if and only if there is a probability measure p (called the Berger measure of T) defined on 
the n-dimensional rectangle R = [0, a\] x • • • x [0,a n ] where = ||Tj|| such that 7^ = j R t^dp{t) := 
J^tf 1 ■ ■ ■ tfrdp(t), for all k e t\. 

Consider now a subnormal 1-variable weighted shift W a , with Berger measure £, and let h > 1. If we 

let 

M h := \J{e n :n>h} (2.2) 

denote the invariant subspace obtained by removing the first h vectors in the canonical orthonormal basis 
of £ 2 (Z + ), then the Berger measure of W a \M h is ^t h d^(t). 

An important class of subnormal weighted shifts is obtained by considering measures /x with exactly 
two atoms to an d t\. These shifts arise naturally in the Subnormal Completion Problem ([7], [S]) and 
in the theory of truncated moment problems (cf. [6], [9]). For to,t\ £ R+, to < Hi and po,pi > 0, the 
moments of the 2-atomic measure \x := po5t Q + piSt 1 (here 5 P denotes the point-mass probability measure 
with support the singleton {p}) satisfy the 2-step recursive relation 

7„ +2 = <A)7n + v?i7n+i (n > 0); (2.3) 

at the weight level, this can be written as a^_i_i = + (n > 0). The atoms to and ti are the zeros 
of the generating function 

g(t) := t 2 - (pjt - <po. (2.4) 

More generally, any finitely atomic Berger measure corresponds to a recursively generated subnormal 
weighted shift (i.e., one whose moments satisfy an r-step recursive relation); in fact, r = card supp p. 
In the special case of r = 2, the theory of recursively generated weighted shifts makes contact with the 
work of J. Stampfli in [23], in which he proved that given three positive numbers uq < a\ < a 2 , it is 
always possible to find a subnormal weighted shift, denoted Wr^ ai aa y , whose first three weights are 
OiQ, a.% and a 2 . The shift T = Wf aQ (Xi a2 y received special attention in [8], and has a 2-atomic Berger 
measure as above; letting a := Oq, b := a\ and c := a|> we often refer to this shift as the abc shift. We 
will have occasion to use these shifts in Sectional 



3. Statement of the Subnormal Completion Problem 

Definition 3.1. Given m > and a finite family of positive numbers £l m = {(a^, Pk)}\k\<m> we sa U ^hat 
a 2-variable weighted shift T = (Ti,T2) with weight sequences a£ and p£ is a subnormal completion of 
if (i) T is subnormal, and (ii) {ot^iP^) = (ak>Ac) whenever |k| < m. 

Remark 3.2. Note that since a subnormal 2-variable weighted shift is necessarily commuting, £l m in 
Definition 13.11 satisfies the commutativity condition in (|2,ip . When a family of positive numbers has this 
property, we say that it is commutative. 
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Definition 3.3. Given m > and a finite family of positive numbers O m = {(akj Ac)}|k|<mj we sa V that 

(l m +i = {(«k 5 /3k)}|k|<m+i * s an extension of Q m (ak, AO = («k>/?k) whenever |k| < m. T/ie degree 
of Q, m , deg fi m; is m + 1. When m = 1, we say that Oi zs quadratic. For m = 2£ + 1, the moment 
matrix of fi m is 



M(f) EE M(J2 m ) EE M (O r 

Observe that if Q m +i is commutative, then so is Vt m 



'— (T(*J)+(p,ff))°| i +^ m • 

Form odd, M(Vt m+ 2) is an extension of M(0, T 



Notation 3.4. When m = 1, we shall let a := 



00' 



6 := 



/^00> 



a 



10) 



d := 



Oqi and / := 



To be consistent with the commutativity of a 2-variable weighted shifts whose weight sequences satisfy 
(|2.ip . we shall always assume af = be. The moments of fii are 




= b 
be 



720 := 6d 



and the associated moment matrix is 




In this case, solving the SCP consists of finding a probability measure /i supported on such that 
f R2+ y % x 3 dfj,(x,y) = 7^ (i, j > 0, i + j < 2). 

Associated with the measure /i of a subnormal completion is the moment matrix 



M(2)[ M ] := 



/ 7oo 


7oi 


7io 


702 


7n 


720 \ 


7oi 


702 


7n 


703 M 


712 M 


721 [fA 


7io 


7n 


720 


712 M 


721 N 


730 [fA 


702 


703 [fA 


712 [fA 


704 [fA 


713 M 


722 [/i] 


7n 


712 [fA 


721 M 


713 M 


722 M 


731 M 


\ 720 


72i Em] 


730 M 


722 M 


731 M 


740 M / 



(cf. dTUl). 



If we now let 


p : = 


«20; := 




:= an d 




/ 


1 


a 


b 


ac 


6e 


6d \ 








a 


ac 


be 


acp 


beq 


bdr 








b 


be 


bd 


beq 


bdr 


bds 








ac 


acp 


beq 


704 [fA 


713 M 


722 M 








be 


beq 


bdr 


713 [fA 


722 M 


731 M 






V 


bd 


bdr 


bds 


722 [fA 


731 M 


740 Em] / 





Q2, we see that 



The (quartic) moments of [i give rise to an extension of Q\, so that M(2)[/i] = M(Qs). It is thus 
clear that a necessary condition for the existence of a measure n is the positivity of M (O3), which in turn 
implies the positivity of M(Oi). If we now let p := d^o, 9 := "ii> r := "02 an d s := /3 ( 

M(0 2 

The localizing matrices M X {VL^) and M y {Cl^) (cf. [11, Introduction]) are 

fa ac be \ 
Mr (^3) = ( oc acp 6eg | and Mj,(fi3 




(The matrix M^f^) is the compression of M{Vt^) to the first three rows and to the columns indexed by 
monomials containing X, that is, X, X 2 and YX; the matrix M y (Cls) is defined similarly.) Observe that 
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M x (Q,s) = M( 0j i)(l) and M y (Cls) = M( 1>0 )(1) (cf. (jl.ip ). For the existence of a measure n supported in 
M+, it is necessary to have M x (Q,s) > and M y (Cl^) > 0. 

In this paper we prove that starting with the positivity of M(Oi) alone, it is possible to choose new 
weights p, q, r and s to ensure the positivity of M X {Q,^) and My{Vt^). We can do this while simultaneously 
building a positive flat moment matrix extension M{Vt^) of M{Q\). Once we establish the simultaneous 
positivity of M x (0,^) and My(0,s), the existence of a representing measure [i follows from the 

main result in [TTJ . We prove this in Section In Section [6] we give a concrete description of p in terms 
of the initial data a, b, c, d and e and the new weights p, q, r and s. First, we present in Section U] an 
abstract solution to SCP, which uses our new approach, involving localizing matrices and the results in 

mi. 

While the flat extension approach is successful in the case m = 1, it will not lead to a solution of SCP 
in all cases. Indeed, it is possible to build a moment matrix M(2) = M(Q^) admitting a representing 
measure, but with no flat extension M(3) (cf. Section [7] below) . This shows that our approach, while 
very general, will not yield subnormal completions merely by one-step flat extension techniques. In 
many instances, solving SCP will require a finite sequence of rank-increasing extensions followed by a flat 
extension; this is despite the fact that for SCP one looks for a measure with support in the nonnegative 
quarter-plane. As a matter of fact, the "translation of support" technique we use in Section [7J shows 
that solving SCP is equivalent to solving A-TMP, where K is a compact set satisfying K C and 
K n (0, +00) 2 7^ 0. Thus, SCP is a special case of A-TMP, and it is natural to expect that qualitative 
aspects of TMP theory will be appropriately reflected in SCP. 

4. Abstract Solution of SCP 

In this section we will give an abstract solution of Problem 11.31 We first consider the main theorem 
in [Tl]. Although Theorem 1.6] deals with truncated complex moment problems, there is an entirely 
equivalent version for the case of two real variables, which we now state. 

Theorem 4.1. Let V = {p\, . . . ,ppj} C C[x,y] and define fc, by deg pi = 2fc, or deg pt = 2ki — 
1 (1 < i < N). There exists a rank M (n)- atomic representing measure for ^( 2n ) supported in 
K-p := |(x,y) E M. 2 : pn (x, y) > 0, 1 < i < A} if and only if M (n) > and there is some flat exten- 
sion M(n + 1) for which M Pi {n + kj) > (1 < i < A). In this case, the representing measure for 
M (n + 1) is rank M (n)- atomic, supported in K-p, and with precisely rank M{n) — rank M Pi (n + kj) 
atoms in Z := {(x,y) G M? : Pi(x,y) = 0} (1 < i < A). 

With the aid of Theorem 14.11 we can now state and prove a result which gives a sufficient condition 
for the solubility of SCP in two variables. Our version does not completely match the conditions listed 
on Theorem [L2l and we now explain why. In one variable, building a flat moment matrix extension of a 
Hankel matrix entails adding an extra row and an extra column, and checking that the rank is preserved. 
This entails checking the range condition in Theorem ll.2( iii) and ensuring that the new lower right-hand 
corner entry satisfies the requirement in Smul'jan's Lemma |22j : 

Lemma 4.2. (cf. [91 Proposition 2.2]) Consider the 2x2 block matrix D := f ^ q^) ■ Then 

D > A > 0, B = AW for some W, and C > W* AW . 

In two variables, what one adds is not a row and a column but instead a block of rows and a block 
of columns; while it is still possible to preserve the range condition, the new lower right-hand corner is 
not a number but a square matrix, which must necessarily be Hankel for the extension to be a moment 
matrix. One easily finds out that ^-hyponormality (cf. Theorem |1.2( ii)). while necessary, is no longer 
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sufficient to prove the hankelicity of the new lower right-hand block. Thus, our result avoids mention 
of £-hyponormality. Moreover, solving the SCP admits two structurally different cases: m odd and m 
even. In the former case, deg Vt m {= m + 1) is even, so we have enough moments to build the moment 
matrix M{Q. m ). 

The same is not true, however, when m = 2k, since we have moments up to degree 2k + 1, and this 
does not allow us to build a complete moment matrix. In the terminology of Lemma 14.21 we have 
A := M(O m _i), and also the B block (consisting of moments up to degree m + 1), but no C block. Since 
we are seeking a moment matrix M(Q m +i), with moments up to degree 2m + 2, we can certainly require 
that Ran B C Ran A = Ran M(O m _i), but that in itself does not generate the additional moments. 
One could attempt to define the C block as W*AW (where W solves the equation AW = B), but this in 
general does not produce a Hankel block C, as has been observed in |12| . Therefore, it becomes necessary 
to postulate the existence of moments of degree m + 1 that, together with the initial data O m , allows us 
to build a moment matrix, which we will call M(£l m+ i). 

Theorem 4.3. Let Q m := {(«k> Ac) : |k| < m} be an initial set of positive weights satisfying the commu- 
tativity condition Ac+ei^k = «k+£ 2 Ac k S with |k+£j| < m (i = 1,2)), and let m := 2 [y] + 1; 
thus m = m if m is odd and m = m + 1 if m is even. Assume that M{Vtfn) > 0, and that VLfh admits 
a commutative extension such that the moment matrix M(£lfh+2) is a flat (i.e., rank-preserving) 

extension of M(£lm), with M x {VLfh+2) > and My(Clm+2) > 0. Then there exists a rank M {Vtfh)- 
atomic representing measure \i supported in R^, with precisely rank M(Qrh) — rank M x ((lm+2) atoms 
in {0} x M + (resp. rank M(£lm) — rank M y ({lm+2) atoms in IR+ x {0}). The measure /j, is the Berger 
measure of a subnormal completion Vt^ of Sl m , provided at least one atom of /i lies inside the positive 
quadrant in M?. 

Proof. In the case at hand, the polynomials pi are p\{x,y) := x and P2{x,y) ■= y; thus, k\ = ki = 1. It 
follows that K v = M.\ and that M P1 (n + ki) = M x (n + 1) and M p2 (n + k 2 ) = M y (n + 1). Our result 
now follows from a straightforward application of Theorem 14.11 □ 



Despite its simplicity, Theorem 14.31 is quite useful, as we will see in the next section. We conclude this 
section by showing how the additional moments required in case m is even are sometimes determined by 
M(O m _x). 

Example 4.4. Let m = 2 and assume that A := M{Q\) > and det^4 = 0. Then there exist moments 
7i,j (* + j = 4) sucn that M(f^3) > is a flat extension of A. The case when rank A = 1 is easily 
disposed of, so without loss of generality we focus on the case Y = al + bX in the column space of A. 
We are assuming that A > 0, M x (tl 3 ) > 0, M y (n 3 ) > and Ran B C Ran A. (Observe that M X (Q 3 ) 
and M y (Q 3 ) include moments up to degree 3, so building them requires no new moments.) The equation 
det A = uniquely determines 702, from which we obtain at once the weight 



a oo/^oo a io 2 a oo/j?oo a oi + PpQaQ! 
«oo(«?o -«oo) 
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Since Ran B C Ran A, each column in B must be a linear combination of the columns 1 and X, and 
straightforward calculations using Mathematica yield unique values for 020 , "n and 002- Concretely, 



a 



20 — 9 / 9 9 \ 

«Io(«oo - «0 



'017 

2 _ ago«f aoi - aooapi - aoo« 2 «o2 + 2 «oo«oiQo2 ~ ^01^02 

11 ™2 f„,2 ^,2 1 



«oil«oo - a 0lJ 



,2 



Q Oo(^00 Q! 10 /^00 a 01 + a 00/^02 a 10^02) 

^oKo-O 

With this information at our disposal, it is now straightforward to check that the C block, defined as 

C := mW*AW (where AW = B) is Hankel. Thus, M(Q%) := f ^* ^ ^ is a moment matrix extension 

of A, and moreover rank M(Qs) = rank A = 2. It is now clear that SCP admits a solution in this 
particular case. 

One might wish to extend the above reasoning to the case rank A = 3, as follows. Let W := A _1 B and 
let C := W* AW . It is well known that C is in general not Hankel, and that one can make it Hankel by 

adding a rank-one positive matrix P. Thus, M(Q-^) := p* q ^ p J ^ s a positive moment matrix, 

and rank M(Q^) = 4. The solution of the Quartic Moment Problem [12] now says that there exists a 
flat extension M(VL§) of M{Vt^). Unfortunately, we can't tell whether the support of the representing 
measure for M($7s) is contained in the first quadrant in M?. This would require verifying that the 
localizing matrices M X {Q,^) and M y {n^) are positive. If we knew that they are flat extensions of M X {Q,^) 
and M y (Q^), resp., then of course we would be done. This fact is false in general, but it might be true 
in the context of SCP; however, we have not been able to prove it for SCP. 

5. Localizing Matrices as Flat Extension Builders 

We now specialize to the case m = 1 in two variables, and show that the condition M(f2i) > is 
sufficient for the existence of a subnormal completion. 

Theorem 5.1. Let Q\ be a quadratic, commutative, initial set of positive weights, and assume M(£li) > 0. 
Then there always exists a quartic commutative extension of Qi such that M(i}^) is a flat extension 
of M{Q,\), and M X {Q,^) > and M y (Cl^) > 0. As a consequence, £l\ admits a subnormal completion 
T- 

Proof. Since m = 1, we have 1 = 1. By Theorem 14.31 we first need to show that six new weights, 
«20) P20, $11, «02 and /3o2 can be chosen in such a way that M x (Qs) > and M y {p,^) > 0. Once we 
prove this, we shall employ techniques from truncated moment problems to establish the existence of a 
flat extension M(Cl^) of M(Ox). We will then appeal to the main result in [11]; the existence of a flat 
extension will readily imply the existence of a representing measure fi for M(l), and the positivity of the 
localizing matrices M x {2) and M y {2) means that supp [i C IR^. Thus, /j, will be the Berger measure of 
a subnormal 2- variable weighted shift Tq^, which will be the desired subnormal completion of fix- 
We now build M(2). To simplify the calculations, we let 



«oo b 
af d 

«oi / 



A 2 
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(The family Qi is shown in Figure [TJ) 



Figure 1, The initial family of weights Oi 



Thus, 



1 a b 

M(l) = | a ac be \ . (5.1) 
b be bd 



Since M(l) > 0, it follows that det ( ^ ^ ) > 0, i.e., 



By the commutativity of Q\, we have 



and therefore 



acd > be 2 . (5.2) 
af = be, (5.3) 
cd > ef. (5.4) 



A straightforward calculation shows that 

det M(l) = acbd - b 2 e 2 - a 2 bd + 2ab 2 e - b 2 ac 

and that 

det M(l) > cd- ef > 0; (5.5) 

for, if cd — e/ = then the rank of the 2x2 lower right-hand corner of M(l) is 1, and then M(l) cannot 
be invertible. Inspection of (I5.4p reveals that we must have c > e or d > /. Without loss of generality, 
we shall assume that c > e. We also assume that a < c, since otherwise a trivial solution exists. (In 
fact, if a = c in ()5.1|) . the positivity of M(l) implies that a = e and b = f < d; when 6 = d (resp. b < d), 
the point mass <5( a ,6) is the Berger measure of the subnormal completion (resp. (1 — ^)#( a ,o) + ^(a,d))- 
Thus, in what follows we shall always assume c > e and a < c. 

To build M(2) = M(^), we first need six new weights (the quadratic weights), namely ct2o, P20, 
an, $u, ctQ2 and /3o2- Since the extension ^3 will also be commutative, two of these weights will be 
expressible in terms of other weights. We thus denote 020 by y/p, an by ^/q, &02 by ^Jr, and /?02 by ^fs 
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(/?20 and Pn can be written in terms of the other four new weights 

/ 1 



M(2) 



a 
b 

ac 
be 
\ bd 



a 
ac 
be 
acp 
beq 
bdr 



b 

be 

bd 

beq 

bdr 

bds 



ac 
acp 



be 
beq 
bdr 



bd 
bdr 
bds 



It follows that 
\ 



(5.6) 



(with the lower right-hand 3x3 corner yet undetermined) and 

be 



M x {2) 




beq 
bdr 



and M y (2) 




Now, since the zero-th row of a subnormal completion of f2i will be a subnormal completion of the zero-th 
row of f^i, which is given by the weights a < c, we let p := c. By one of the main results in [18], having 
«io = «20 immediately implies that a\\ = y/c, that is, q := c. Thus, 



MJ2) 




By Choleski's Algorithm [2], M x {2) > if and only if bdr > so that we need r > . Thus, provided 



we take r > the positivity of M x (2) is guaranteed. It remains to show that we can choose s in such 



a way that s>d and M y {2) = M y (2)(s) 



Case 1: 

also let s := 



d. 



M y {2) 



It follows at once that rank M y (2 



M y (2)(s) > 0. We consider two cases, 
c. By (|5.4p we have d > f, so we can take r := c and guarantee that M x {2) > 0. We 
We then have 

b be bd 
be be 2 bed 
bd bed bd 2 

1, and therefore M y (2) > (and of course s > d). 
Case 2: e < c. We define r by this extremal value, i.e., r := ^t. This immediately implies that 

011 := VVi an d by propagation, /3u := \ff (all j > 2) in any subnormal completion. The resulting 
weight diagram is shown in Figure [2j 

It remains to define s, in such a way that s > d and M y (2) > 0. Since 

be 
bee 
bef 



M y (2) = M y (2)(s) 




and the 2x2 upper left-hand corner of M y (2) is invertible, we see that M y (2) > if and only if 
det M y (2)(s) > 0. Since det M y (2)(s) is linear in s, we pick for s the unique value that makes 
det My(2)(s) = 0. A straightforward calculation shows that 

_ a 2 cd 2 - 2abde 2 + b 2 e 3 
a 2 d(c — e) 



We then have 



d 



e(ad — beY 



> 0. 



a 2 d{c — e) 

Thus, this particular choice of s guarantees both s > d and M y (2) > 0. 



n 









V d 






Vd 


77 










Vb 


77 


77 


fa 







Figure 2. The family f2i augmented with the inclusion of the quadratic weights 



To complete the proof, we need to define the 3x3 lower right-hand corner of Af(2), and then show 
that M(2) is a flat extension of M(l), and therefore M(2) > 0. We consider the following two cases. 

(i) rank M(l) = 2: Without loss of generality, we may assume that a < c, so that the columns 1 and 
X of M(l) are linearly independent. The column Y must then be a linear combination of 1 and X, 
and that allows us to define YX and Y 2 in M(2). Moreover, since the zero-th row of is given by 
the weights -fa, 7c, -7c, • • • , whose Berger measure is £ x = (1 — ~)8q + ~5 C (and thus supported in the 
two-point set {0, c}), it is natural to let X 2 := cX in the column space of M(2). With these definitions, 
one easily verifies that the truncations to the first three rows of X 2 , YX and Y 2 agree with the 3x3 
upper right-hand corner of the matrix M(2) in (|5.6I) , It is clear that the matrix M(2) thus defined is 
positive semi-definite, but one needs to verify that M(2) is a moment matrix. This amounts to checking 
that the (4, 6) and (5, 5) entries are equal. Now, a straightforward calculation shows that in the column 
space of M(l) we have 

b(c — e) f — b , 
Y = -i '- • 1 + J - X, 5.7 

c — a c — a 

so that 



f-b 
1 + X, YX 

c — a 



M(2) 
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' Y 2 X 2\ 



[Y 2 ,cX) 



c (Y, YX) = c 



b(c - e) 



c— ^-be + c— bee 

c — a c — a 

bee— = beef. 
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On the other hand, using (pTTj) we define YX := ^r^-X + ^X 2 , so that 

M(2) 55 = (YX,YX) 

b(c — e) f — b o 
-i -X + X 2 , YX 

c — a c — a 

= b -^bce + l^(cX,YX) 
c — a c — a 

= — -6ce + c— bee 

c — a c — a 

= bee— = beef. 

c — a 

It follows that Af(2)46 = M(2)s5, as desired. In this case, the representing measure is supported in the 
two-point set {(0, yo), (c, y c )}, where 

He ~ e) 

yo ■= (5.8) 

c — a 

and 

b(c-e) f-b cf -be 

Vc ■= 1 c = = /. (5.9) 

c — a c — a c — a 

(iii) rank Mil) = 3: We let B denote the upper right-hand corner of M(2), that is, 

ac be bd \ / ac be bd 
B := | acp beq bdr I = I ac 2 bee bdr 
beq bdr bds J \ bee bdr bds 

We also let C denote the lower right-hand corner of M(2). Since we want rank M(2) = rank M(l) = 3, 
we must define C := B T M(iy 1 B. A gain, we need to verify that Af(2)4g = ^(2)55, i.e., C13 = C22- A 
straightforward calculation shows that 

C13 = bedr. 

When c > e, we have r = and another calculation shows that 

b 2 ce 2 



Q 



22 



a 

,2, 



it is then immediate that C13 = C<i2- When c = e, we have r = c, and in this case C13 = C22 = be d, as 
desired. 

The proof of the Theorem is now complete. □ 

6. Description of the Representing Measure 

In this section we provide a concrete description of the Berger measure for the subnormal completion in 
Theorem l5.11 We have already observed that when rank M{1) = 1, the representing measure is /i = <5( a ,&)- 
When rank Mil) = 2 (and the columns 1 and X linearly independent), there is a 2-atomic representing 
measure, with atoms (0, yo) and (c,y c ) given by (|5.8p and (|5.9|) : thus, \x = P(o,y )^(o,y ) + P{c,y c )$(c,y c )- 
To find the densities p(o, yo ) and pi c ,y a ), we use the first two moments: J dp, = P(n yVo ) + P(c,y a ) = 1 
and J s dp, = cpr c yc \ = a. It follows that the densities are P(o,y ) = 1 — t an( ^ P{c,y c ) = §• Thus, 

V = I 1 ~ f)<W) + c 6 (c,Vc)- 

We now focus on the case rank Mil) = 3. Since Mil) is invertible, the last three columns of the 
flat extension M(2) can be written in terms of the first three columns; that is, the columns labeled X 2 , 
YX and Y 2 are linear combinations of 1, X and Y. Each of these column relations is associated with 
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a quadratic polynomial in x and y, whose zero sets give rise to the so-called algebraic variety of |13j ; 
concretely, V(&s) := C\ p (x Y)=o deg p<2 ^>(p), where Z(p) denotes the zero set of p. In our case, the three 
column relations are 

X 2 = cX 
YX = fX 



The associated zero sets are 



^9 be(f — d)^ r cd — ef^ r 
Y 2 = — Y- ±X H -Y. 

a(c — e) c — e 

{(x,y) : x = or x = c} 

{(x,y) : x = or y = /} 

2 be(f-d) cd — ef 
{{x,y) ■■ y = —, r aH y\- 

a(c — e) c — e 

Let z := c ^~g and observe that z > by (I5.5p . The algebraic variety of ^3 is then Vitlz) = 
{(0, 0), (0, z), (c, /)} and these are the three atoms of the unique representing measure for M{2). To 
find the densities, we use the first three moments, 700, 701 and 7io : 

P(0,0) + P(0,z) + P(c,f) = 1 

P(cf)C = a 

P(0,z)Z + P{c,f)f = b. 

a 

c 

2 



We obtain 



P(c,f) 



l f , a b(c-e 
P(o,z) = ~{b-~f) 



z c c(cd — ef) 

P(0,0) = 1 " P(0, z) -P(c J) = aK J_ ef) det M(l). 

Thus, the representing measure is 

det M(l) b(c-e) 2 a 

" = aftCcd-e/)^) + c(d-e/)^» + (<U) 
Direct calculation shows that J s 2 d/i(s,t) = ^c 2 = ac, J st d/j,(s,t) = -cf = af = be, and 

o , , , b(c — e) 2 9 a „ 9 b(c — e) 2 ,cd — ef s9 6ef 

i 2 ^ S) i = ^ 2 + -/ 2 = ^ ^ -f + - 

c(cd — ej) c c(cd — ej) c — e c 

_ 6 M -e/) + te/ =H 

c c 

so that correctly interpolates fii. 

Recall now that the marginal measures and z/*' associated to a Borel measure v on the Cartesian 
product 1x7 are given by v x (E) := v(E x F) and u Y '(F) := v(X x F), for £ and F Borel sets. In the 
specific case of the measure /U in (|6.ip . observe that the marginal measures \i and /tr are (1 — -)Sq + -8 C 

an< ^ ab(cd-^ef) ^° + c(cd-ef) ^ z c^/' res P ec tively. While /U X is always 2-atomic, /x y is 3-atomic if and only 
if £ 7^ /. When is 3-atomic, its moments (which are also the moments of an associated unilateral 
weighted shift W v ) satisfy the recursive relation 7„+2 = —fzj n + (/ + z)7 n+ i (all n > 1), with 70 = 1, 
7X = b and 72 = bd. It is easy to see that the restriction of W v to the invariant subspace Mi defined in 
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(|2.2p has Berger measure -j-i dfjX(t) = P{0 '* )Z 5 z (t) + p< -°'^ whose recursive coefficients are —zf and 

z + /, respectively. 

On the other hand, it is indeed possible to have z = f, which occurs precisely when d = f. In that 
case, the three atoms of n are (0,0), (0,/) and (c, /), and the unilateral weighted shift associated with 
fx Y is W(h : d,d,— )- The reader will note that the location of these atoms can also be predicted by Theorem 
14.31 once we observe that rank M x (2) = 1 and rank M y (2) = 2. 



7. Flat Extensions May Not Exist 

We now present an example of a set Q3 for which the associated moment matrix M(2) admits a 
representing measure, but such that M(2) has no flat extension M(3). Thus, while Theorem 14.31 provides 
a general sufficient condition for solving SCP, not all SCP will fit that framework, and their associated 
moment matrices M(fim) will require a sequence of moment matrix extensions M(£l r n+2), • • • , -^(^m+2fc)) 
with M(£lfh+2k) admitting a flat extension M(0^ +2 (fe + i)). 

The example is motivated by the construction in |12^ Examples 1.13 and 5.6], and also by |12t Propo- 
sition 1.12], which states that a TMP and its image under a degree-one transformation of the base space 
are equivalent as moment problems. In particular, the qualitative aspects of TMP are preserved under 
degree-one transformations; our idea is therefore to "translate" [12, Example 1.13] three units to the right 
and four units up, so that the support of the 6-atomic representing measure in |121 Example 1.13] will 
land in the positive quadrant. (We note that to produce a valid representing measure for SCP, it suffices 
to have all atoms in the nonnegative quadrant, and at least one atom in the positive quadrant.) To 
effectuate the above mentioned translation, we recall the definition of the Riesz functional L 7 associated 
to a TMP. The linear functional L 7 acts on polynomials by L(y l x 3 ) := jij. Given the moments 7^', 
one can translate the TMP by h units in the horizontal direction and k units in the vertical direction by 
letting 7y = L^u 3 ) := L^((v + 4)*(u + 3) J ). The associated moments of degree 4 are: 



Too 



= 1 



701 = 1 7io = 1 

702 = 2 711 = 720 = 3 

703 = 4 712 = 721 = 730 
704- ™~-n 



700 = 1 

701 = 4 

702 = 17 

703 = 76 



7io 
711 
712 



= 5 
19 
77 



720 = 27 

= 4 712 = 721 = 730 = 9 703 = 76 712 = 77 721 = 97 730 = 157 

= 9 713 = 722 = 731 = 740 = 28 704 = 354 713 = 331 722 = 371 731 = 535 740 = 972 



For example, 



721 



LMv + A) 2 {u + 3)) = L 7 ((w 2 + 8v + 16)(« + 3)) 



L 7 (v 2 n + 8vu + 3v 2 + 16u + 2Av + 48) 
+ 8- + 3- 3 + 16 -1 + 24 -1 + 48-1 



721 + 8711 + 3720 + I6701 + 24710 + 48700 
97. 



With the new moments at hand, we form the matrix M(2). The corresponding weights are: 



«03 
«02 

«oo 



V535 

VTFf 
3^ 

V5 

2 



«12 

an 

OT0 



/37T 



97 
77 
10 
17 
2 



(7.1) 



«21 
"20 



/331 



V77 
2v{l9 
17 



«30 



17 

38 
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Example 7.1. Let Q3 be given by (|7.ip and let M(2) = M(2)(J23) the its associated moment matrix, 
with entries built from the data 7^'. Let M(3) be a positive semi-definite, recursively generated, moment 
matrix extension of M(2). Then rank M(3) > rank M(2). As a consequence, M(2) admits no flat 
extension M(3). For, consider a moment matrix extension 



M(3) := 



1 


4 


5 


17 


19 


27 


76 


77 


97 


157 \ 


4 


17 


19 


76 


77 


97 


354 


331 


371 


535 


5 


19 


27 


77 


97 


157 


331 


371 


535 


972 


17 


76 


77 


354 


331 


371 


705 


714 


723 


732 


19 


77 


97 


331 


371 


535 


714 


723 


732 


741 


27 


97 


157 


371 


535 


972 


723 


732 


741 


750 


76 


354 


331 


705 


714 


723 


706 


715 


724 


733 


77 


331 


371 


714 


723 


732 


715 


724 


733 


742 


97 


371 


535 


723 


732 


741 


724 


733 


742 


751 


157 


535 


972 


732 


741 


750 


733 


742 


751 


706 / 



where the moments of degree 5 and 6 are new. A direct computation shows that rank M(2) = 5, and 
that (X — 3)(Y — 4) = 0, that is, YX = 4X + 3Y — 12. In any positive semi-definite, recursively generated, 
extension M(3) this column relation would still be valid, and it would also give rise to two new column 
relations, namely YX 2 = AX 2 + 3YX - 12X and Y 2 X = 4YX + 3Y 2 -12Y. These three identities 
lead at once to the values 714 = 1497, 723 = 1513, 732 = 1925, 741 = 3172, 715 = 243 + 4705, 724 = 6555, 
733 = 7375, 742 = 10796, and 751 = 1024 + 3750- Now, since the compression of M(2) to the rows and 
columns indexed by 1, X, Y, X 2 and Y 2 is invertible, we can find coefficients A\, A x , Ay, A X 2 and A Y 2 
such that 

A X \\\ B + A X [X] B + A Y [Y] B + A X 2[X 2 ] B + A Y 2[Y 2 ] B = [X 3 ] B , (7.2) 

where [-] B denotes the compression of a column to B := {1, X, Y, X 2 , Y 2 }. A calculation using Mathe- 
matica [21] reveals that A x = -25513 + 157 05 , A x = 13587 - 8705, Ay = 1, A X 2 = -1692 + 705 and 
A Y 2 = 0. If M(3) were a flat extension of M(2), an identity similar to (|7.2p should hold for the last row 
in M(3), that is, 

A 1 [l] {x , } +A X [X] {X3} +A Y [Y] {X3} + A x2 [X 2 ] {x : i} + Ay 2 [Y 2 ] {x3} = [X 3 ] {x - iy . (7.3) 

Using Mathematica again, it is easy to check that ^4i[l]{x3} + Axp^jj^a} + ^y[^]{x 3 } + A X 2 [X 2 ]^ x sy + 
j4y2 [F 2 ]{x3} = 7376, while P^ 3 ]{x 3 } = 733 = 7375. It follows that M(3) cannot be a flat extension of 
M(2). 

Remark 7.2. The SCP in Example 17.11 does admit a solution, and the subnormal completion has a 
6-atomic Berger measure. We see this after we observe that the positive semi-definite moment matrix 
extension M(3), while not a flat extension of M(2), does admit a flat extension M(4). Rather than 
showing the details here, we refer the reader to \12\ Proposition 5.5 and Example 5.6]; the representing 
measure constructed there must be translated three units to the right and four units up to give rise to 
the Berger measure that solves SCP in Example 17. 11 
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